Slow regions percolate near glass transition 
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A nano-second scale m situ probe reveals that a bulk linear polymer undergoes a sharp phase 
transition as a function of the degree of conversion, as it nears the glass transition. The scaling 
behaviour is in the same universality class as percolation. The exponents 7 and (5 are found to be 
1.7 ± .1 and 0.41 ± 0.01 in agreement with the best percolation results in three dimensions. 

PACS numbers: 64.60. Ak, 64.70.P, 82.35.-x,83.80.Sk 
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I. INTRODUCTION 

Polymeric melts, or linear polymers in the bulk, ex- 
hibit a rich array of different behaviours, as a function 
of the temperature, specific volume (per chain), chain 
length (molecular weight), and timescale of observation, 
generally without sharp boundaries between the different 
regimes. Linear polymers in the bulk range from solid, 
glassy or leathery materials at low temperatures to rub- 
bery, and eventually viscous fluids at higher temperatures 
or polymer to solvent volume ratios [0^1). The viscosity 
of bulk linear polymers typically increases with increasing 
chain length and/or degree of conversion as polymeriza- 
tion proceeds, and diverges near the glass transition 
The glass transition temperature is a function of the aver- 
age chain length, for small to moderate length chains, but 
ceases to depend on the molecular weight for relatively 
long chains. The rheological or viscoelastic properties 
can thus be studied experimentally either as a function 
of the temperature, or as a function of the reaction time 
during the conversion process, as we have done here. 

In this paper we would like to report that the extremely 
short-time response of bulk linear polymers probed by a 
fluorescent aromatic molecule displays a sharp transition 
as a function of the reaction time, which is a (nonlinear) 
measure of the degree of conversion. Immediately after 
the transition point, the melt is found to be glassy or 
rubbery for temperatures respectively below and above 
the glass transition temperature quoted for commercially 
available bulk PMMA g, namely 106°C. 

We find that this transition is characterised by perco- 
lation critical exponents. In an unpublished paper H], we 
interpreted, erroneously, the percolation exponents found 
here as signalling the onset of the entanglement percola- 
tion transition. We now believe, for reasons to be dis- 
cussed in this paper, that the fluorescence signal cannot, 
in fact, detect the formation of entanglement clusters. 

In order to understand the physical nature of the pro- 
cesses underlying this percolation transition, one must 
follow the reaction kinetics, compare results with ex- 
periments directly measuring rheological properties in 
the course of the polymerization reaction, and also criti- 
cally examine the steady state fluorescence technique em- 
ployed. The "rheo-kinetics" of bulk polimerization is an 



extremely active and growing fleld of research, where dif- 
ferent considerations based on free- volume, entanglement 
(or reptation), gel-effect and vitriflcation compete with 
and complement each other. |^ We will interpret our re- 
sults in the light of recent developments in the theory of 
the glass transition. 0, ^, ||, |l^ 

In section 2, we outline our experimental method, in 
section 3 we discuss the various mechanisms which come 
into play as the reaction proceeds. In section 4, we ana- 
lyze the results in terms of percolation theory and extract 
the scaling exponents. In section 5, we present a discus- 
sion of the results. 



II. THE EXPERIMENT 

Our experimental method, the in situ monitoring of 
free-radical polymerization by a steady state fluorescence 
technique, is the same as the one we used to determine 
the critical exponents for the gel fraction (/3) and the 
average cluster size (7) in a previous paper where 
we also found percolation exponents. 

Fluorescence measurements yield direct, in situ in- 
formation regarding the connectivity and available free 
volume in polymer melts and gels. The use of flu- 
orescent probes to study the different environments in 
which the probe molecule flnds itself, including exper- 
iments on polymerization chemical gel forma- 
tion ijl^, p ^ , swelling |^ and slow release of large 
molecules ||20||is an extremely well established method. 
Such measurements have also been applied to the study 
of glass formation, especially by the group of Ediger [pi[. 



A. Experimental procedure 

Pyrene was used as the fluorescence probe to 
detect the free-radical polymerization process of 
poly(methymethacrylate) (PMMA), a linear polymer. 
The characteristic time for the direct relaxation by fluo- 
rescent emission of the pyrene molecules is ~ 10~^ sec- 
onds, a time scale much smaller than those involved in 
reptation, and of about the same order as Rouse dynam- 
ics II, H, |2|, |23, l2|, |6l. Thus, extremely fast "snap- 
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FIG. 1: The sharp increase in the fluorescence intensity ob- 
served as a function of time for different temperatures. The 
maximum intensity is normalized to 100. 



shots" of the effective network can be obtained as the 
polymerization reaction proceeds, without mechanically 
disturbing it, as one would have to do with conventional 
viscosity measurements. 

Our experimental setup is the same as in ref. ex- 
cept for the absence of a crosslinker, so that only lin- 
ear polymers are obtained. The free radical polymeriza- 
tion of MMA was performed in the bulk in the presence 
of 2,2'-azobisisobutyronitrile (AIBN) as an initiator, at 
60°, 70°, 75°, 85° and 125° C. The glass transition tem- 
perature of bulk PMMA is 105°C for the average chain 
lengths attained here. The pyrene molecule was excited 
at 345 nm and the variation in the fluorescence emission 
intensity, /, was monitored with the time-drive mode of 
the spectrometer, with a frequency of up to 10 points 
per second, staying at the 395 nm peak of the pyrene 
emission spectrum. 



B. Results 

In Fig. 1, we present our results for the normalized 
pyrene fluorescence intensities as a function of the reac- 
tion times, for different temperatures. The results for 
125°C are very similar; however due to the much greater 
mobility and hence higher reaction rates at the higher 
temperature, the transition occurs already at 300 sec- 
onds, and therefore is off-scale in this figure. 

The end product is rubbery at 125°C and glassy oth- 
erwise. In any case, at the end of the polymerization 
reaction the bulk polymer does not flow, but retains its 



shape. The addition of solvent causes the polymer to 
swell, as a gel would, up to volume ratios of around 1:1. 
This indicates that the sample has definitely undergone 
a transition to a gel phase, due to the physical entangle- 
ments of the linear chains. Larger quantities of solvent 
dissolves the top layers, until, at a polymer to solvent 
ratio of about 1 : 2.5 ± 0.2 the whole sample becomes 
liquid. Since we started off with a bulk polymer, this ra- 
tio corresponds to the volume ratio at which the chains 
disentangle. This is in the same ballpark although a bit 
bigger than the conversion estimated from the value of / 
at the inflection point in Fig. 1, frompe = le/Ima.^ — 0.3 
(see [p7|)and the percolation threshold for a cubic lat- 
tice in three dimensions {pc = 0.31) [Q. For a cubic 
lattice, Kantor and Hassold |^ find an "entanglement 
threshold" very close to the ordinary percolation thresh- 
old, pe~pc-2x 10-'^. (Also see H |oHl|. It should 
be noted, however, that these authors consider clusters 
including branch points as well as those that cannot be 
pulled apart because of topological constraints.) 

These considerations seem to point in the direction of 
a gelation (entanglement) transition |3|, |3^ 
at the onset of the sharp rise, ending, however in a glassy 
phase, within a very small portion (duration wise) of the 
whole polymerization process. The glass transition tem- 
perature is swept up as the polymerization proceeds, un- 
til it crosses the temperature of the sample. p5| Even in 
the case of the run with the temperature held at 125°C, 
one is within about 20 degrees of the glass transition tem- 
perature, and one has to take into account the precursor 
effects that may be present. Qualitatively very similar 
behaviour observed many tens of degrees above the glass 
transition temperature |15| , remains, however, a puzzle. 



III. REACTION KINETICS AND RHEOLOGY 

A close examination of the reaction kinetics of the 
free-radical polymerization process ^ re- 

veals that in the initial stages (the flat, vanishingly small 
parts in Fig.l, preceding the rapid rise) the release of the 
free-radical which initiates chain formation, and the self- 
or cross-termination of the chains, balance each other 
out. Thus the chains which are formed are all more 
or less of the same length; i.e., as conversion proceeds 
at a fairly constant rate, the lengths of the chains that 
are formed are not growing on the average. As more of 
the monomers are converted, we reach a point where the 
fluorescence intensity starts picking up. This coincides 
with the rise in viscosity of the medium, with a concomi- 
tant rise in the rate of reaction, with the trapping of the 
free radicals]^, termed the "gel effect." (Some 

authors refer only to the point where the reaction rate 
reaches a maximum, as the "gel effect" |3^). Somewhat 
counterintuitively, the reaction rate grows as the reac- 
tion becomes diffusion limited, because end to end- or 
self-termination of the chains become increasingly rare 
and the length of the growing chains can become very 
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large as a result. We will argue in this section that the 
fluorescent signal grows because of the trapping of the 
pyrene probes as well as the free radicals, in regions of 
relatively high density of the polymeric chains, where the 
cage effect comes into play. ||2^, |4^, ^ 

The curves in Fig. 1 are quite typical of polymeriza- 
tion reactions with or without ||ll|, ^ ||, |o[ H 
crosslinkers, and display a very sharp rise, with an inflec- 
tion point, where the time derivative of the curve passes 
through a maximum, as shown in Fig. 2. (The dura- 
tion and slope of the initial region of low conversion may 
depend on the temperature and weight fraction of the 
crosslinker or initiator.) 

In an experiment involving a free-radical co- 
polymerization with cross linkers. Okay et al ]39| ] have 
monitored the viscosity of the medium using dilatomet- 
ric techniques, and found that the gel point is located 
roughly at the onset of the gradual rise in the pyrene in- 
tensity curve, and deflnitely before the inflection point is 
reached. In the present experiment, since no crosslinkers 
are used, the increase in the viscosity may be ascribed to 
the physical entanglement of the polymer chains pO[ , giv- 
ing rise to a purely geometrical gel phase (note that the 
unconverted monomers are easy to displace and provide 
the "sol" phase). The PMMA chains that are formed 
interact very weakly with each other, so that inter-chain 
attraction as in the formation of thermoreversible physi- 
cal gels Q is less important here. 

While chemical gelation in the presence of chemical 
cross-linkers exhibits a sharp transition which is mod- 
elled by percolation theory 46, 4^, Q, the formation 
of physical gels by bulk linear polymers is less well under- 
stood. It may be conjectured on the basis of scaling ar- 
guments and renormalization group calculations ||2|, ^ 
that this gel point is also controlled by percolation ex- 
ponents, although it arises purely from entanglement ef- 
fects and there are no (chemically bonded) branch points 
in the spanning network thus formed. Nevertheless, the 
above dilatometric flndings on a similar system (albeit 
having crosslinkers) suggest that the pyrene fluorescence 
is probing the post-gel region and not the vicinity of the 
gel point. Further evidence for this is provided by a crit- 
ical examination of the mechanism for the unquenching 
of the fluorescence itself. 

The steady state fluorescence technique relies on the 
fact that the relaxation of excited aromatic molecules 
via indirect non-radiative transfer of energy to the high- 
frequency vibrational modes of the small solvent parti- 
cles or monomers is forbidden |50t] with the increased 
viscosity of the medium. The rates of intra-molecular or 
inter-molecular non-radiative transitions are determined 
by the availability of appropriately spaced levels either 
in the molecule itself or in the surrounding molecules. 
These levels are altered with the degree of conversion of 
the polymer melt, due to the local dynamical constraints 
arising from a reduction in the available free volume, and 
the consequent slowing down of segmental motions, lead- 
ing to vitrification. |2^, 52, 5^, ^ |5^ In this context, 
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FIG. 2: The typical critical peak exhibited by the derivative 
of the fluorescence curve, at 75°C, in the vicinity of the crit- 
ical point. The maximum is at tm, the half width at half 
maximum, fi/2, and the critical point, tc are indicated. 



free volume will mean the volume occupied by uncon- 
verted MM A molecules. 

The connectivity may also contribute indirectly to re- 
tarding the local dynamics, in the present case through 
entanglement effects (since no chemical cross linkers are 
present). But the crossover between Rouse to reptation 
dynamics takes place above the glass transition (in tem- 
perature, or free volume). Moreover, it is generally be- 
lieved that the cage effect Q leading to this slowing down 
and vitrification occurs at much smaller spatial scales (in 
fact of the order of the polymeric bonds) js^ than those 
of Rouse or reptation ||2^, |2^ dynamics. 

We can make a rough estimate of Afc, the average 
length of the chains at the inflectionpoint, from the value 
of the conversion fraction pe — 0.3 |27]| , assuming that the 
number of chains is approximately the same as the num- 
ber of initiators. The molecular weights of the MMA and 
AIBN molecules are 100 and 164, respectively, and the 
experiments were performed with samples consisting of 
1ml MMA (density 0.96 gr/ml) and 0.0026 gr AIBN. We 
find Mc ^ 180, which is of the same order of magnitude 
as the chain length at which the chain dynamics crosses 
over from Rouse-like to reptation behaviour p5| , 
as estimated by Kreer et al. from simulations. This 
is an order of magnitude larger than an alternative esti- 
M| which can be obtained by comparing the 
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average end to end distance to the diameter of the rep- 
tation "tube." 

On the other hand, the time scales involved in the cage 
effect are of the order of pico seconds, while in Rouse 
or reptation motion they range from nano to microsec- 
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onds ||2^. Therefore it would seem that those degrees of 
freedom which would be frozen out with respect to the 
pyrene probe with a time scale of the order of nanosec- 
onds, would be the latter. Thus it is difficult to unequiv- 
ocally identify the contributions to the unqucnching of 
the pyrene fluorescence, as arising strictly from the con- 
nectivity (entanglements) or from the formation of den- 
sity inhomogeneities that are the precursors of glass for- 
mation, merely from considerations of time and length 
scales. 

Further information is provided by the ESR spectra 
in the course of the polymerization (free radical copoly- 
merization of MMA and ethylene glycol dimethylacry- 
late - EGDMA) These measurements reveal that 

the radical concentration exhibits a sharp peak right at 
the inflection point of the polymerizations curves, which 
are exactly similar to those appearing in Fig.l. Even 
more crucial is the fact that, in the region between the 
inflection point and where the curves flatten off, the ESR 
spectra switch from having 13 lines to 9 lines; with the 9 
line spectra being contributed by radicals flnding them- 
selves in solid, i.e., glassy environments. |4|, |5|] 
A simulation ||57|| based on a random walk percolation 
model |4^, p9[ compares resonably well with the ex- 
periment below the percolation point, but fails to do so 
beyond the gel p oint, occuring at about weight 30% poly- 
merization. p7| 

These considerations lead us to conjecture that the 
critical scaling behaviour which we observe in the vicin- 
ity of the inflection point of Fig. 1 , and which is discussed 
in the next section, arises from the percolation of dense 
regions where the cage effect can be observed, and vit- 
rification sets in. This also raises the possibility in Ref. 
[ prf , the critical scaling that is observed is due to the 
percolation of slow regions in this sense and may lead us 
to reasses our interpretation of the latter results. 



IV. PERCOLATION OF SLOW REGIONS 

We would like to cast the quasi-static properties of the 
polymer melt in the language of percolation, to interpret 
our results. 

Note that free monomers are relatively easy-to-displace 
objects which do not constrain the motion of the large 
chains. By a slow region, we shall mean a region of suf- 
ficiently high density, such that parts of several chains, 
or even of the same chain, get in each other's way and 
inhibit each other's motion on a short time scale. 

We take the volume fraction occupied by the total 
number of monomers incoporated into the chains as the 
"occupation probability," p, of the sites of a three dimen- 
sional lattice. We will take a cluster to be a set of occu- 
pied points on the lattice, which are nearest neighbor to 
at least one other member of the set. Therefore such 
nearest neighbor occupied sites may be considered as be- 
longing to the same cluster, regardles of whether they are 
chemically connected (i.e., belong to the same chain) or 



not. A spanning cluster of this description will be called 
a percolating cluster. Wc then argue that the pyrene 
molecules contributing to /, the fluorescent intensity, are 
precisely those which are trapped within interstitial re- 
gions of these clusters. 

If M is the average molecular weight of the chains. 



p = MCvo/V = vo/v 



(1) 



where C is the total number of chains, vq ~ V/N is 
the volume occupied by a single monomer, and N is the 
total number of monomers. The average length M of the 
chains obeys the differential equation 



dMjt) 
dt 



k{N - CM) 



(2) 



over a sufficiently short time interval where the time de- 
pendence of the reaction rate k due to the change 
in the viscosity, and of C (due to the gradual dissolu- 
tion and termination of the initiators) can be neglected. 
Thus, for relatively short time intervals, we get a linear 
growth law. 



M{t) - M{to) = [M^ - Mito)] kC{t - to) 



(3) 



where Moo = N/C and t^ is an arbitrary starting point. 
Defining tc as the time at which the entanglement perco- 
lation threshold Pe is reached, and using Eq.(|l|), we have, 
for sufficiently small |i — <c|, 



|t-tc| 



1 \P-Pe 
kC 1 - Pe 



OC |p - Pe 



(4) 



Since below tc, there is no percolating cluster, the total 
normalized fluorescent intensity will be proportional to 
the average cluster size S. For t > tc, most of the pyrene 
molecules are trapped in the macroscopic network of slow 
regions, and I then measures the fraction Poo of the 
monomers that belong to the macroscopic cluster. 
For a system with linear size L, the scaling forms for the 
quantities S and Poo around the percolation threshold, 
together with (0) yields. 
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t < tc 
t > tc 



(5) 



Here, (3 and 7 are the critical exponents for the strength 
of the infinite cluster and the average cluster size, and 
df = (/3 + 7)/t/ is the fractal dimension of the span- 
ning cluster, with 1/ being the correlation length expo- 
nent. Notice that we need not subtract the value of 
I{tc) from I{t) in (^) for t > tc since we are assuming that 
once the threshold has been crossed, the unquenched flu- 
orescence intensity is being contributed essentially by the 
monomers trapped in the incipient infinite cluster. 

The time derivative of the intensity I is plotted in Fig. 
2, and looks like a typical critical peak, with rounding 
due to "finite size" effects. The fits to the double log- 
arithmic plots of the fluorescence intensity v.s. |t — ic| 



TABLE I: Experimentally determined values of /3 and 7. 



T(°C) 


60 


70 


75 


80 


85 


125 


7 


0.40 


1.8 


1.6 




1.71 


1.41 


/3 


0.45 


0.42 


0.41 


0.41 


0.41 


0.41 



for t > tc and t < tc, are shown in Fig. 3. The criti- 
cal point tc and the exponents (i and 7 are determined 
by simuhaneously fitting the intensity data to the be- 
haviour in Eq. (|^) , in such a way as to obtain the greatest 
range in \t — tc\ over which one finds scahng, both above 
and below the critical point. |]lT| In this way it is pos- 
sible to remove the ambiguity from the position of the 
critical point and greatly enchance the accuracy of the 
critical exponents. We see that the critical region is con- 
fined to the interval 3 < \t — tc\ < 30. The values of 
the exponents are given in Table I. The agreement with 
the known values (/3 = 0.41,7 = 1.80) of the percola- 
tion exponents in three dimensions, is very good, with 
/3 = 0.41 ± 0.01 and 7 = 1.7 ± 0.1. The error bars are 
estimated from the values reported in the table. The 
numbers for 60° C differ markedly from the rest and we 
have excluded them from our estimate of the exponents. 
At such low temperatures the mobility of the MMA are 
reduced jl^ to the point where the reaction ends before 
all the monomers are exhausted (the absolute intensity 
is lower at this temperature). The remaining monomers 
act as if a solvent were added to the bulk, increasing the 
lower length scale of the percolating network [|ll[, and 
distorting the scaling behaviour. On the other hand, at 
125°C, the sample is boiling before the polymerisation 
has had time to progress, and the curve below the tran- 
sition point is rather noisy, so that we find 7 = 1.4; we 
recover (3 = 0.41 above the transition point. 

The critical point tc obtained from the best scaling fits 
is smaller than tm marking the maximum of the dl/dt 
curve. The quantity c = {tm — tc)/ti/2 where ti/2 is the 
half- width of the dl /dt curve at half maximum, is found 
to be 



tr 



0.248 ±0.005 



1/2 



(6) 



over the whole range of experiments, and for different 
temperatures. This value is the same as that found 
from the measurements performed near the percola- 
tion transition found from fluorescence measurements on 
crosslinked gels [pT| . 

Finite size scaling predicts that this ratio is in fact 
a universal constant for percolation. From Eq. we see 
that dl /dt must obey the finite size scaling relations 



dl 
'dt 



7^(7+1)/!' 



-t)Li/-] 
tc)L'/''] 



t < tc 
t > tc 



(7) 



with the scaling functions (f>s{x) ^ x^'' and (j)p{x) ~ x^ 
for a; ^ 1, and ^ const, for a; <C 1. Let us define the 
universal scahng function <!>' by dl/dt ^ - 




FIG. 3: Best double logarithmic fits to the curve for 75° C 
shown in Fig. 2, yielding 7 = 1.60 (lower curve) and /3 — 0.41 
(upper curve). Every tenth data point is plotted for greater 
clarity. 



tc)L^^'^] ■ Then <&' clearly corresponds to the curve de- 
picted in Fig. 2, both above and below tc, and c must 
be given in terms of its moments, z„ = J z"'^'{z)dz, as 
c = zi/{z2 - zlf/^ - 0(1). The fact that we find the 
same numerical value for c as that determined from the 
cross-linking gelation experiment |pT| is further indica- 
tion that the transition we observe here is in the same 
universality class as percolation. 

To understand this universality, we note that a dense 
region will look like a branching point under a coarse 
graining by some factor 6, (while a bigger "blob" will now 
look like a tight entanglement region) . If the percolating 
entanglement cluster is self-similar at length scales i ^ 
where ^ is the correlation length diverging like \p—pe\~'^, 
near the percolation threshold, this transformation will 
simply renormalize |p — Pe| by a factor b^/", but other- 
wise leave the scaling behaviour unaffected. Below the 
percolation threshold, the finite clusters themselves are 
fractal for £ ^ and the same argument applies. 



V. DISCUSSION 

We have observed that steady state fiuorescence exper- 
iments in a polymeric melt reveal a sharp transition in 
the narrow crossover region from low conversion to high 
conversion and vitrification. The transition is charac- 
terised by percolation exponents. On the basis of argu- 
ments presented in section 3, we believe that the critical 
scaling behaviour we observe in the fluorescent signal as 
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a function of time is not due to the physical gelation of 
the sample due to entanglements, but to the percolation 
of dense regions where the cage effect constrains the local 
segmental motions of the linear chains and the radicals 
as well as the pyrene probe. 

It is relevant to note in this context the dynamical scal- 
ing hypothesis which asserts that, as the glass tran- 
sition is approached from above, long range dynamical 
spatial correlations build up between particles of fragile 
glass formers. Dynamical clusters are formed out of those 
objects which have cooperatively rearranged, using the 
same free volume to accomplish the move, g |, |lO| In the 
present context, this idea acquires particular significance, 
because with the onset of the "gel effect," and the trap- 
ping of the free radicals, further polymerization must pro- 
ceed in exactly the same way as the motion of a particle in 
a vitrifying medium, leaving behind, moreover, a perma- 
nent "trail" consisting of a long chain that is formed as a 
result. The chain can only grow if unattached monomers 
(which act like free volume with respect to monomers al- 



ready incorporated into the chains) can exchange places 
with segments of neighboring chains trapping the radical 
at the head of the growing chain. Once a monomer is 
attached to the growing chain, a new monomer has now 
to exchange places with the surrounding chain segments, 
and will in turn be eaten up by the radical, etc. This pro- 
cess then leaves behind a growing chain confined within 
dense regions, and interconnecting these dense regions. 
The sudden rise in the rate of polymerization at the on- 
set of the "gel effect" feeds on itself in nonlinear fashion, 
and suddenly gives rise to the percolation of dense, slow 
regions in the melt. The entrapment of the pyrene probe 
in these slow regions gives rise to the critical behaviour 
of the fluoresent intensity. 
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